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. TOMOGRAPHY: 
THREE-DIMENSIONAL IMAGE RECONSTRUCTION 



1. THE PROBLEM* 



Suppose we want to mapr'out the interior of an opague 
object. To be specific suppose we want to determine the 

.interior of part of a human body — an arm^ say. ,An x-ray 
photograph reveals the two-dimensional proj ection ^of the 
arm on a viewing screen . .'This projection indicates points* 
of high and low mass density. Jfowever, given a point a£ 
high density* we still do not know where* in the interior 
of the arm the point is located— it could be anywhere on 
the line perpendicular to the screen and °]jence parallel to 
the x-ray beam. The problem is an important orfe. A 
le6ion in\1fhe brain will register on the Viewing screen; 
and it 'is important to know just where a-long the line 
parallel to the x-ray beam the^ lesion is located. So the 
problem is in constructing a jthree-dimensional map from 
two- dimensional pictures. (Actually, the phrase mass 

' density -is used in a suggestive sense *here; it is only one 
factor affecting x-ray impenetrability.) » 1 
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Figure 1. Arrows indicate direction of x-ray beam on a 
circular object. The dark area inside the object, corresponds 
to a region of hlgff mass density. AH three have the same 
projection on, the viewing screen, but the locations of the 
high mas£ density are quite different. » It mjght even be uni- 
fdrnfly spread along the l4n£ parallel* to the x-ray beam as in 



'As we shall see, if x-ray photographs* are taken at * 
many different angles, then the information from all of 
them can be synthesized into the desired three-dimensional, 
map. The purpose of this unit is to show how this is*. 
done / 1 




Figure 2. Two x-ray project ions indicate_the point of 
high mass density. ^ 

, Although the goal is to construct. a three-dimensional 
map, the problem is actually a two-dimensional one. If x, 
y,z-axes are situated so that the x-ray beams are perpen- 
dicular to the z-ax'is, it is enough to construct a map of 
each cross section parallel to. the xy plane., Thus we will 
he concerned solely with a two-dimensional o.bject in the 
xy *pl an 6 with incident x-rays* in this plane. 

. ■ • i ' -. 

2. NOTATION 



Let' f(x,y) be the mass density of the object and let 
<b+j be the angle between the x-axis and jthe incident beam 
of x-rays. Draw u,v-axes as in Figure 3, that-is, the 




Figure 3- Arrows i*n upper left'Nndicate the direction 
of x-rays. 



v-axis points in the direction of the' x-ray source while 
'the u-axis is 'perpendicular, to the v-axis with the u,v-axes 
in tihe »same orientation as the x*?y-axes. Finally; we 
Refine, the ray-sum R(u,<|>) a's the total mass registered u 
unfits from the v-axis wh'en <J>+j is the angle Cf the x-ray 
beam.* That is* R(u q ,4>q) is- the total mass on line t in 
Figure 4(b). * « . ^ 





Figure 4(a)c Graph of R(u,<J>) .for fixed <J>. (b) Line I 



is Uq from origin. 



'Now We assume that all ray-sums R(u,<J>) (for all u 
and all 4>) are known, ,(In actuality we would know the 
function R(u,<J>) for a finite number of angles <fr as we 
direct beams from different directiqns^at the object; and 
•R(u,4>0 would be approximated by interpolation ( for all 
other values of <f>.) Notice; that ft(u,*+ir) ■ R(-u,<j>) (Why?); 
so It is only necessary to find R(u,^) for 0 < <t> < it. 

The definition for fixed u Q and <£ n of R(u n ,4> n )-is « 



.related to the .mass density f(x,y) by 
V*0 



R(u n ,<fr n ) = total, mass along* i in Figure 4(b)' ' 



f (x,y)ds 



where s is arc-length along £.< 
course, just u 



Uq = constant. 



(1) 



u = x COS 
V 



"0 



+ y sin 



p 0 



x *sin 



+ y cos 



The equation of I is, of 
But trigononieVry yields 



x = ii cos 4> Q - v sin 4> Q 
y = u sin 4> Q + v cos <{>q. 



r 



: ' i H J 

Exercise \, Derive these equations: Let P be a point with coordi- 



nates 
Show 



(x,y) in the x,y coordinate system and (u,v) in the u,v system, 
that (x,y) and (u,v) are related as above.- 
v 

▲ 

P 




So 

* G'2) R(u 0 ,* 0 ) 



f(x,y)ds 



f (u Q cos <|> 0 -v<sin 4> 0 *u Q sin <j) Q +v cos 4> Q )dv. 



(Note "that if the object has £iiiite extent [so that 
f(x,y) = 0 for (x,y) noi in the object], then the above 
^ integral actually has fini-te limi/ts,) The proolem is to 
solve this' integral equation for f(x,"y) given the ray-sums 
fc(u,4>)/ In other words, in a real situation the ray-sums 

woulci be known, but the actual mass density 'would -be of 

* & 
interest. * 



3. BACK PROJECTION: AN APPROXIMATION 



Lefts' first obtain -an approximate solution -to the ^ 
problem tff reconstructing the unknown mass density f(x,y) 
from the known ray-sums R(u,4>). This is the technique of 
back projection. , Suppos*e in addition to the ray-sums we 
also know that the object is contained in the disk of 
radius L about the origin. Now R(u,<|>) is the integral. .of 
f along % ; and the length of t (or rather the length of 
th e int ersection Of A with the disk of radius L) is 
2/L 2: u 2 (see Figure 5). Suppose that (x Q ,y 0 ) is a point 



A. 



J — 



* ft 2 

S Figure 5. I has length 2/U -u . 

tf • " ' . * 

on Under the assumption Chat all the mass on I is 

uniformly distributed 'along the line segnfent I, then the 

mass density almig.il wpuld be + 

(3) 6(x y )•=' . ' * 

f . 2/L 2 -u* . . 

where u =* x Q cos <p + y Q sin 4> as in Exercise 2. Finally 
we approximate f( x 0> v o) a ^ the average of all ^(Xq,Yq) as 
<J> paries from 0 to 2tt. 

\ ' 

« Exercise 2 . Show, that the equation of line i through (xAy A ) at . 
- & % s* o o • \ 

angle <J> as in Figures is ' % 

* * • . • 

• *• - y ^.-(tan^)" 1 x + [y Q + (tan <f>)V x Q ] . t> 

. Show that the di statue from this line to tfie origin is 

/ u = x o cos y o sln 



/ 

Using Exercise 2 we see .that the back a .,p,roj ection esti- 
mate of f is • • 



t * R( x n cos ♦ *< yn* n < f l > < f l ) 
■ —J j—^ j-d* . 

, a 2 If, 2 - (x 0 cosV + y 0 sin<f.) 2 l 

Example : Suppose the object is actually the disk of radius 
L of uniform density f(x,y) *v constant c. Then we would 
* observe the ray-sums (from Equation' (2) ) 




\ 



Thus for aevery u 



• ' ■ R K») I c for |u| < L 
and back projection yields^ approximation 



' F <V>V. = 27 



2tt 



cd<p - c • 



t° T . ^o ,y oV insid £ * he di sk of radius L.^ And Yhis agrees 
with the *exact value of f(Xg f y 0 J. 



Exercise 3 . Suppose the object is contained in the disk of radius L \ 
ahout the origin with density f(x,y)*> c(x 2 + y 2 )i-c a positive 4 * con- 
stant. .Find the ray-sums R(u,<J>) and the back projection F(x,y).*»How 

does?'F compare with' f? • a 

• * i * 

^ ExeVcise k . Suppose the. object 'is. known to be contained, in the di9k 
of radius L about the origin,, but is actually the disk of radfos 
L Q < L of constant density c. 
■» c 



AbscPssa represents 




J 



r =* »/x 2 +y < " = distance from origin.- 



*What are the ray-sums^ R(u,<{>) and the back projection approximation 
F(x,y)? ' (Leave your ^answer as a* integral.) Draw a rough sketch of 
F as a function of the distance r fromrtne originT* 



In general the back projection approximation over- 
estimates regions of low mass density and underestimates 
regions of high mass density. This is due t j.o the fact 
that projection fo,r each fixed $ of R(u/<J>) uniformly back 
through the object will overestimate regions of low 
density and^ underestimate regions of high density. 

4. FOURIER TRANSFORM NOTATION , 

> Let f(x) be a function of the real variable x,.$j 

-c* < x^ < ». We use this convention for the Fourier.' Trans- 
it 

form f of f: j • ' * 



£(«) • 



f(x)e' 27Tia,x dx, 



Under suitable conditions the Fourier Inversion Theorem 
holds ■ 



' f(x) - 



^ r ^ 2rtiu>x, 
f(a>je do>. 



(Suitable ^conditions are, for' example, that 



f 2 jA * <, 

|f| < (», f isy||||tinuous and is piecewise differenti 



able.J) These results can easily be extended to functions 

f(x,y) of two variables. Thus we set 

I , 1 .... / 

f(o>,o) =^ 



/ 2 2, » 
Exercise 5 > Let, f (x>y) = e~** x Y ; where'a Is a pos ? 1 1 ve ^pstant , 



U$e*contaur integration to show* that 

2 



' f(« f o)'* i«cp J [- (a) 2 + o 2 )], 



Show that the Fourier I nyefs ion . Theorem holds for this f . (Note: 
the second part of this problem requires^no new'integrations,) 



5, THE , EXACT SOLUTION 

We now show 1 how the ray- sums R(u,<J>) completely deter- 
mine f(w,o), tjie Fourier Transform of the 'mass density. 
The trick is* to change the variables of integration in the 
integral defining f above. Switch to the u,v coordinates 
defined in Equations (1).. Orient w,o axes — the domain 
variables^ for( f — so that <J> = arctan £ . Then 



- 1 7 -U 

cos <J> = sec <J> = (1 + tan <J>) ^ 

= (1 ♦ a 2 /* 2 )'* " 



/(i) +0 



s-fn <J> -= 



• a 



n — 2 



wx+oy = w(u cos (J) - v sin <J>) + o(usii\4> + v cos <J>) 



And under suitable conditions the Fourier Inversion Theorem 
holds * 



-Also dudv » dxdy since '(u',v) and (x,y) are each Cartesian 



f 



ERLC 



T(x,y)' 



4 

12 



PS 



dudo . 



'<.<,' coordinate systems:.. So 



f (u,o) 



(5) 



'f(x,y)e- 2,Ti l ux+0 y)dxdy' 



f(x,y)dv e*p 



^-2iri(/ 



a) +o )u du, 



v. i 



JL O 



"where x - u cos $ - v sin <j> 
x y = u sin <J> + v cos . 

But the inner integral is R(u,<J>) by Equation^) . Thus 



f(a),o) = 



exp 



•27Ti( 



/> 2 +a 2 )u 



du. 



However, phh§ integrS^ is by (definition the-.Fourier Trans- 
form of R(u,$) with respect to .the u variable at /u> 2 +o 2 
and with $ fixed at arctan (cr/w) . Irf other words, f(w,a) 
is obtained by fixing > dt arctan (oAo) , taking the* 
Fourier Transform of R(u,tf with respect to the single" 
variable u, and evaluating this Fourier Transform at 
/I — I 

/w +o . For a function g(x,y) of two variables we denote 
the Fourier Transform with respect to the first variable 
by g^Casy) : ( * 



g x (w,y) = 
Summarizing yields 



g(x,y)e- 2 ^ iwx dx. 



6. THE BASIC THEOREM 



(6) 



f(o>,a) = R x 



/(i "~ ------ 1 



a> +o , arctan — 



The Theorem shows us how, in 'theory, to reconstruct the 
mass density f(x,y) from the ray-sums R(u,<J>). Namely,' from 
R(u,<J>) calculate the Fourier - Transform with respect to u, . 
RjCx, <)>)•. Then f(o>,o) can be calculated by the Theorem from 
which f(x,y) can be found by Fourier Inversion. Of course 
in practice R(u,<f>) is much too complicated to find its 
Fourier Transform without a cpmputer approximation of the 
integral defining R 1 (t,4>) and similarly for the step from 
f to f. Still, the Theorem^ndicates the general technique. 



Example : Suppose the ray-sums are found to be 

2 ' . 



R(u,4>) = e 



•au 



for- positive constant a. (Noting that R(u,<J>) is indepen- 
dent of <J>, we expect f(x,y) to be a function of the 

2 - 2 

distance x +y from the origin only.) *Now 



-au -27T1TU, 

e e du 



exp' 



f 2 , 

TT 2 

T T 



either using tables or by *a contour integration as in 
Exercise 5 * So 



,f(a),o^) = R^(/w +o , arctan (o/w) 



i ex P 



— ((^ +0 ) 



Thus 



f(x,y) = 



c e v 2irHo)x+ay) , ,5 
f(w,a)e +v '^dwda 



2tTia)x 



dio 
! 



exp 



TT 2 2 A ~ . 

— o + 2TTiay 



do 



-^4 



/a -ax 

TT 



'a -ay 

TT C 



% exp 



-a(x 2 + y 2 ) 



1 r \ 
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where the second, to the last equality $>liows/by a'ccjmtour 
'integration as in Exercise. 5/ 1 y ' j 



/ 

-fe , 



Exercised . , Su PP° se *(u,<f>Hs found to be e\' U ' independently !of <J>. 
Find Rjd,^), f(w,o) A , set up* the integral for f(x,y). ; * 



1'. ClkcULAR SYMMETRY - " { 

, - i 

i< - i " i 

Ma.ny corollaries follow from tfte basic theorem <{>f 

Equation (6). We indicate on^ whic^h generalizes Exejcise 

6 and the Example before it. k : /' , % [" 



Suppose R(u,£) ' exhibits, some/sort of periodic' ' 



behavior *in <J>. That is, suppose /the x-ray photograpns 
begin 'to repeat themselves ars we/move the 5 angle of the 
incident x : ray fceajn. (Of cour*s^ such a repetition ta^jways 
occurs by moving the beam -by 2nL We have smaller periods 



of repetition in mind here.) then it is natural to use 
jit • i ~ 

Polar *rather thin Cartesian . Coordinates . Thus let • 



V 

x = r- cos 



r sin 



as always for Polar Coordinates. And in the (<p,o)* plane 

a) = t cos B 
o = t sin B 

denote the Polar representation. .That is, t ■ /w 2 +cr 2 , 
3 a 'aTctan' (o/u>) . Now 



C7) f(x,y) 



f(a»,o)e m ^ x+ °yid^ * 
V 

R^/^7 2 ", arctan (o/ Je 27rilwx+ °y^da)do 
A ° -11 



Rir 



2n 



■R 1 (t >B )e 2llirtC0S ^"^tdBdt 



where the second equality follows from the Basic Theorem 
arrd the last by using Polar Coordinates: the fact -that 

• wx+oy =.rt cos 6 cos 6 + rt sin 6 sin 6 

X = rY cos (6-6) v . \ 

.Note that tdgdt is the area element -in polar coordinates. 
And b>5 Exercise 7 below we conclude that 



(8) 



2tt 



R 1 (t,3 + 6)e 2TrirtC0S ^td6dt. 



Exer yL se 7 ' Let 9 b e a real valued function of a real variable which 
is periodic of period 2 : 

9(&+27T) = g (B) 

fox all ,3- Let c be a fixed number. *By considering the graph of g 
show that ^* i 



27T+C 



g(3)d6 = 



2tt 



-g(6)dB. 



Noting that both cos 6 and Rj(t,3+9) are periodic of period 2tt 

(so that their product is as well) deri ve\£quation (8) from the last 

expression in Equation (7). ~ 



So far we/have been -completely general. To understand 
the above equation mor^let's take the case R(u,<J>) = R(u) 
independently of <J>. Then R^t,^) = R ][ (t) is independent 
of <J> (Why?). Sfo * * 



(9) 



f(r,9) 



2tt 



fit* 



1 7 
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Noting tnat the integral is independeivC of 8, we "have just 
proved the obvious fact: If the ray-sums R(u,$) are 
independent of the angle of the incident x-raysj then f 
-must be circularly symmetric y i,. e . , independent of 8. 

Exerciss 8 . let X be a real constant. The value df 
tin 

\\ cos s, 
e ds 



Is denoted 2ttJ 0 (A). (J q is called the Bessel Function of order zero.) 
Evaluate this integral as an infinite sum using contour integration. 
(Hint: Convert to a contour integral over the counterclockwise upjt 
circle centered at the origin using z =,e ,S . Show that the integral 



x 

which is (2TTi) Residue of exp|i j 
ate this residue as an i nf in i te %er ies . ) 



at z = 0. Now eval u- 



. Using Exercise 8 and its solution we see that i£ 
fc(u,<J>) is' independent of <J>, then (from Equation (9)| 

(10) m f(r,8)' = f (r) 

r 2ir 



Ri ft) 



0 



e /irirt cos S d6 



0 



tdt 



R^Ct) •.2iTJ 0 (^TTrt)tdt t 

0 ' 



where JqU) 'is ihe Bessel Function of order zero, 

I 2j 



V x,.. FJ^ltf. 1 . 



. (J!)' 
■ X 
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(See the Answer to Exercise .8.) This formula exhibits 
f = f(r)^as what is called the Hankel Transform of order 
zero of R^t)*. Hankel Transforms occur quit/* naturally in 
two-dimensional Fourier Analysis problems. ' * ** 



Exercise 9 . Suppose the ray-sums are found to be - independent of <f> 
with. 

'• • r . < " 



R(u,<f>) = R(u) = 



1 + (2ttu) 2 



By a contour integration show that R^(t)*= e'^L Show that 

r 

_t ^t 2j + , dt » (2j + 1)! 



e 
0 



and that f(r,6) = 2tt(1 + (irr) 2 )" 372 . .(Note: You will need the fact 
that w * * 

an.d the Binomial Theorem in the special case 

1 1 

(. + x)- ,/2 = I xJ 

* * j=o 1 J > 

where by def ini tion 

J-l/2j „ (-1/2)(-1/2-l) „ . (-1/2- \ *\) ) 



Although^this' formula exhibiting f (r) as the Hankel Trans- 
form of RjtO solves the qase when R(u,<£) = R(u) is 
S-ndependentrof *4>, thje actual integration would normally 
/require a computer approximation. However-, the Hankel * 

/Transform of some* common functions arellisted in mathematical 

'tables. • " V* 



14 
» 



/ 



. f 8, THE CIRCULARLY. ASYMMETRIC CASE 



A similar analysis can be carried out when R(u,<j>) 
doe% deptend on (J).' The solutioii exhibits *f(r,e) as combi- 
nations* of higher order Hankel Transforms. The solution 
is elegarit, but we relegate it to • 



Exercise 10 . NowRj (t ,<J>) — the Fourier Transform of fi(u,<J>) with 
respect, to jthe u variable— is for each fixed t a periodic function of 
• $ of period 2tt. This is so since"the ray-sums R(u,<J>) repeat as <f> 
changes by 2tt. Thus Rj(t,<J>) can be 'expanded , as a # Fourier Series for 



each fixed t: 



in<J> 



where by the Four ier Series formula 



C n«> " Iff 



2tt 



1 



a. Show that c (tj Is the Four ier Transform of the nth Fourier 
n — 

Coefficient of R(u,<J>) 



<n< U > " H 



2ti 



R(u,4.)e '"'d*. 



b. Plug the new expression for ( t »4>) into Equation (8) to show 

that . " * - • 

t • ' 

* .-co 

int) 



t 

where 



c (t)d (rt)tdt 

n T n 



4) 



27T 



exp [i ns + 27tirt-cos s] ds. 



By definition J (A) = (27ri n )~ £\U/2tt) is caHed the Bessel 
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c. Now fix r and consider f(r,6) as a function ^>f, 6. Conclude from 

10a and 10b above, that the nth Fourier Coefficient of f(r,0) 
' when expanded in a Fourier Series, as a function of 8 is 



c (t)j (27rrt)tdt 
n n 

- * ' 

which is i -times the nth Hankel Transform of c^. 



- HISTORICAL AND y gLIOGRAPHICAL NOTE 

Analytic .techniques in tomography were first used in 
radioastronomy . Antennae were unable to focus on the 
individual points of the solar surface, but were able to 
focus on thin striplike segments crossing the sun. Thus, 
in making a map of the microwave radiation density of the 
sun, astronontexs had to reconstruct it. from "strip sums." 
In medical radiology x-ray- CAT (Computer assisted tomog- 
raphy) . scanners were introduced in 1972. These machines*, 
take x"-ray pictures at about "one .degree intervals between 
0° and 180°. While CAT scanners are expensive (of the 
order of Ijalf a million dollaVs* per machine), they afford 
a revolutionary improvement in medical diagnosis. 

f%r an interesting nontechnical source of information 
see "Computerized Tomography 1 * by W. Swindell and H. Barrett 
in Physics Today (December 1977). For a deeper mathe- 
matical treatment see "Principles of Computer Assisted 
Tomography 11 by R*. Brooks and G. DiChiro in Physics in^ 
'Medicine and Biology, Volume 21 (1976, pp. 689-732). v For 
an applied, very readable account of Fourier Series, 
Fourier Transforms and Hankel Transforms see The Fourier 
Transform' and its Applications by R.M. Bracewell (McGraw- 
Hml, 1965). An informative, nonmathematical article on 
a similar, method of image reconstruction is "Ultrasound 
in Medical Diagnosis 11 by fieyey and P. Wells found in 



Scientific American <May 1978, pp.- 98-112), In this 
technique-, sound rather than i-rays are used to avoid 
organ damage, 

* 10. ANSWERS AND SUGGESTIONS TO EXERCISES 



3. RW*) - y /l 2 -u 2 (L 2 +2u 2 ) from Equation (2), 

c 2 2 2 
F ^ x 0 ,y 0^ = 3^ L +x 0 +y 0^ from Ec ? uation 



f?. 2 

As functions of r = /x +y the graphs o y f f and F are 



C-2LV3 - 




<u R(u,4>) 



.o, 



L Q < |u| < L, 



So by Equation (2) 



F( V'V "H 



A 7 - 



where u ■ x Q cos (J) + y Q sin <J> and the range of integration is 

over all 0 < 4) < 27T so that lul < L A . Note that 

— — — 0 " 

cL 



F(x Q ,y 0 ) - 0 



n — f 

for »x 0 +y 0 = L. A rough sketch of F as a function of the dis 
tance from the origin is 



>0 "» 



0 

ERIC 




f(w,a) = 



• At 



e -^i( U x + ay) e -a(x 2 + y 2 ) dxdy 



-27Tt(ox -ax . , 
e e dx< 



^27Tiay -ay , 
e 'e ' dy. 



But the one-dimensional Fourier Transform of e" ax ,s 
A f* 2 2] £ 

/J W " X Ur m tables or b Y contour integratic 
Now wl\^ f(w,o) = j exp£- ll (co 2 +o 2 )j, 



f(x,y) • 



tt rir 



8 TT 2 /a 



exp 



-= — (x.+y ) 
tt / a 



1 



i e -a(x 2 + y 2 ) 



by the above integrat ion* wi th — replacing a. So 

a a a 

f(-x,-yj = f(x,y) = f (x,y)—which is the content of\the Four 
Inversion Theorem. » 



-R T (a,*) = 2(l+iTT 2 T 2 )"" 1 .' ' 
Thus 

f(u),o? = ^i+^V+o 2 )!* 1 . 

On the unit ciro^e 



Z*= e is , cos s -i(z+ i), ||- ds. 



Now 



a'(z-i) _ I W , f (a:- 1 )" < • 
i=o j; k=0 . k! 



The coefficient of in the product series ie 



00 J i • 00 2j 
Y il il - V a J 

j=o J - J - j=0 (j!)*. , 



Thus the coefficient of z 



.1 f. A 



T I 
j=0 



[i>) 2j V_J__i y Hii fAl2j 



2rri times this is the answer by the Residue Theorem*. 

. 9.; With R(t) = e ' the formulas afcove the statement of th 
problem imply < 



f (r,6) .= 2tt 



•t . \ H£ (Trr) 2j t 2j + i d ; 
j-o ■ 



2* I (*r) 2j • (2j + l)! 

j-0 (ji) 2 

2tt £ ?j (-l)J' (2j + l)(TTr) 2 j 

*to!ti(V]h)V)v;. 

*d?r-i(T]K) J 
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1. GENETIC DISEASE 

v 

1.1 Introduction 

Genetics, the science of ' heredity, is barely one 
hundred years old. Although people have long been aware 
that parents pass many of their own traits on to their 
children, it is only recently that scientists have under- 
stood exactly how and why this happens. In the process 
of learning the mechanisms of heredity, they have dis- 
covered a large number of gertetic* diseases which, like 

other physical characteristics, a<£ passed down from 

■ 

one generation to the next. 

The gene is the unit of heredity.: A complex protein 
molecule called deoxyribonucleic acid, or DNA, determines 
the color of the eyes, the texture of the hair, and mo'st 
of the other traits that distinguish one person from an- 
other. In addition, ge^nes control a multitude of bio- 
chemical reactions that take place inside cells. When 
they fail to do this correctly, the result can be a serious 
disruption of ^normal metabolism. 

1.2 Cystic Fibrosis and Huntington^ ..Disease 

Cystic fibrosis, the most common inherited killer 
of children*, is caused by just such an inborn error of 
metabolism. This error causes cells to produce a thick 

c 

mucous secretion that clogs the passageways of many vital 
organs. Although treatment for 'cystic fibrosis has im- 
proved in recent years, few of its victims live to be 
adults. If a*couple has one child with cystic fibrosis, 
the chances that each of their other children will be 
born with the disease are one out of four. 

Another genetic dfsorder , Huntington's disease, does 
"not ~makTTtT ^appearance untilthe mid<TTe adultf years 
*At that time, its victims exhibit od,d pastures, involuntary 
motions, and bizarre mental changes. Eventually they 
succumb to the disease.. Many victims of Huntington's dis- * 

1 
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ease have children before they discover that they are ill 
When this happens, the children have one chance in two of 
developing the disease themselves. 



2. MECHANISMS OF INHERITA NCE 
* — i 

/■ 

2.1 Recessive and Dominant Diseases 

To understand why the probability associated with t 
Huntington *s disease is higher than the one associated 
with cystic fibrosis, it is necessary to understand the* 
different ways in which these two diseases are inherited.. 
Gene|, come in pairs of two. In reproduction these pairs 
split, and the child receives one member of each pair, 
chosen at random, from each of his parents. 

• > *-s. 

Each gene has a variety of -different forms, called 
alleles. Certain genetic diseases are caused by a single 

-abnormal- allele. Some of these diseases occur only in 
people who have- two copies of this abnormal allele, one 
from eacft parent. These are called recessiVe diseases. 
Cystic fibrosis is a disease of • this ' type . Others occur 
in people who have only one copy of the abnormal allele. 

'If either parent passes down this abnormal allele, 'the 
child will be afflicted. Diseases like this are said to 
be ddminant. Huntington's disease is an example of a 
dominant disorder. Once it is known that a hereditary 
disease is either recessive or dominant, simple laws of 
probability can be used to determine the chances that 
parents with particular combinations of genes, called 
genotypes, will transmit the disease to their child'ren. 

2.2 Genotypes and Their Probabilities 

A child with cystic fibrosis has two abnormal al- 
leles-, -one -from- each -parent . Thus, if two healthy parents 

a' 

have^ a child with a, recessive disease. like cystic fibrosi^ 
both*of them must carry the abnormal allele. If a symboli 
zes this abnormal allele and A symbolizes the normal one, 
the parents 1 genotypes must therefore be Aa . People\^who 

. - 2 

u4 



; ^ 

have one normal allele and one abnormal one for 'a reces- 
sive disease are called carriers; they usually s*how no 
signs of the disease, but they can pass the abnormal 

allele on to their children. A chil.d with a recessive 

i 

disease has genotype aa. The diagram in Figure l % shows 
how* parents with genotype Aa can have a child with geno- 
type. aa, 

Allele s ^om Mother , 



Allele 
, from 
Father 











A 


a 


A 


AA 


Aa 


a 


aA 


aa. 






*- 



Figure 1, 



Genotype of 
Child 

The possible genotypes of a child with two 
carrier parents. 



Because the parents transmit each of their genes, 
with equal probability, the four possible outcomes are' 
equally likely. There is therefore one chance kn .four' 
that ar child of these parents will get two normal alleles, 
two chances in four, ^r one in two, .that he will be a 
carrier, and one chan<^ in four that he will be -born with' 
the disease. These probabilities are of little use in 
predicting the outcome of any one partici^Lar birth. How- 
ever, of 1000 children born to parents who are both car- 
riers of cystic fibrosis, approximately* 250 caa .be expected 
to be victims of the disease. Of the remaining 750, about 
500 will be carriers. The. 'chances that the healthy child 
of two carrier parents is himself a carrier are therefore 
5^0 out of 750 or two out of three. 



Because it is a dominant disorder, Huntington's, dis- 
ease will manifest itself in people with only one abnormal 
allele, that is, in people of genotype^ Aa, Since Hunting- 
ton's disease is very rare, its victims are almost always 
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nuirri£d xo people with two normal alleles. As .shown in ^, 
< Figure 2, the ghildren of such a marriage*- will be^f w twfc. 
genotypes, Aa „and AA, • (In this case the mother is "the -v. 
normal parent and the father is the carrier.) 

Allele from Mother. 
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from 
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Father 
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Figure 2, 



tfhe possible genotypes of a cjiild with one 
carrier parent. , 



■ ) 



Since Aa individuals will develop the disease an<f AA indi- 
viduals will t not, the. chances are two out of four, or one 
out of two, that the child of a victim of Huntington's 
disease will eventually develop it himself, 

' 3, GENETIC COUNSELING 



3,1 Jpredicting the Occurrence of Genetic Disease 

Although few hereditary diseases can be»treated ef- 
fectively at present, tests have been * developed to detect 
carriers of many genetic disordois* and! to identify a % num- 
ber of defects in unborn children. By this kind of % 
screening, and by ass.es*sing the liklihood that genetic 
disease will arise in particular families, many hereditary 
diseases can Yiow.be prevente'd. Central .to this effort is 
the genetic counselor, one of whose Jobs it is to calcu- 
late the— risk that- particular prarrenrs wMi trails ml^a^ 

hereditary disease/ To do this, r she relies on two functa- 
mental principles of probability, 

1, Probability has no memory/ jhe genotype of each 
, . o i 4 



- - child of^a particular marriage is "independent 
. *of the genotypes -of any previous children. 

^Thus, jthe chances ' that the child \>f two carriers 
-^ of cystic fibrosis will be afflicted with the 
dis-ea-se are one in four "for every c hjft d re'gard- 
~v les;s °* tne outcomes of previous births. 

two or\*nore events are independent, their , 

... probabilities are multiplied to get the proba- 

i-. "-bility that they will occur in sequence. 
-1 * * " * 
• 3.2 Some Typical Problems , 

Here are some problems of the .sort that v a genetic - 
counselor might be called, upon to solve.' ; 

Both' Mr.' and Mrs: B. had -sisters afflicted with 

cystic fibrosis. What is the probability that they'will 

have a child with the disease? I 

y 

Solution : Because they produced afflicted children, 
both Mr. B's parents and Mrs. B's parents were carriers 
of cystic fibrosis.. Of the normal children of carriers, 
two out of three are carriers themselves. ?herefore, 
% the probability that Mr. B is a carrier is 2/3 and the 
probability that Mrs\ B is a carrier is alsq 2/3. Since 
the probability that two carriers of cystic fibrosis wilL 
produce an afflicted child is 1/4, the probability . that ' 
Mr. B is a carrier , ^Mrs . B is a carrier, and their child 
will be born with the disease is 2/3 x 2/3*x 1/4 or 1/9, 
It is worth noting that -the probability of the complemen- 
tary event, that their child will not be bofcri with the 
disease, is 1 - 1/9 or 8/9. ( # * 

The mother of Mr". L's father die'd of Huntington's 
disease. What are the chances that Mr. L will' get it? 

Solution : pie probability that Mr. L 1 s. grandmother 
passed the abnormal allele on to her son rs 1/2. The 

probability that he in turn passed it on to Mr. L is 
-also 1/2, So the probability that Mr. L will develop 

Huntington's disease. is 1/2 x 1/2 or 1/4. 



Exercise 1, Both Mr. L and Mrs. I had brothers aff licted.with .cystic 
fibrosis. If they have two children/ what is the probability that 
both will be afflicted with the disease? 

Exercise 2 . Mr. and Mrs. P arc* tested and both are jfound to be 
carriers of Tay-Sachs disease, a rare recessive disorder. What is 
the probability that their first child will be bom with the^isease? 
That their, first two children will be healthy? That their first 
three children will be carriers? 

Exercise 3 . Albinism is a recessive disorder characteri^d by^ a 
marked reduction in pigmentation throughout the body. A woman who 
is an albino marries a normal man whose brother is an albino. What 
are the chances that their first qhild will be an albino? If their 
first child is an albino, what are the chances that their second 
' child will be too? 

Exereise 4. Two first cousins marry. If the. sister of their grand- 
father had cystic fibrosis, what are the chances that they will have, 
a child with the ^disease? 

3 . 3 Condi tional Probability 

Example 1. Both Mr. and Mrs. A had sisters who died 
4 of cystic fibrosis. If thSA's have two children, what 
are the chances that neither will be afflicted with the 
disease? , .' 

There are two possibilities to be considered here. 
Unjess Mr. atid Mrs. A are both carriers, their children 
run no risk of getting cystic fibrosis. ■ If, on the- other 
hand, both Mr. A and Mrs. A are carriers, the chances that 
each of , their children will be healthy are three out of 
four. Because both of these conditions must hie considered,, 
this problem is said to involve conditional probability^ 

In the diagram in Figure 3, C represents the event 
that Mr.^and.Mrs. 0 A are both carriers, and C is its com- 
plement;. Since both of their parents were carriers, Mr 4 .. 
A^and Mrs. A each have a 2/3 chance of being carriers r them- 

v On 6 

\ KJ J 



v- < 

selves (see Section ^2), The probability of event C is 
hence 2/3 * 2/3 or 4/9, and the,probabil i ty of C is 1 - 4/9 
or 5/9, * ' , 4 

H stands for the event that Mr. and Mrs. A have two 
healthy children, H for its complement. If C is true, the 
probability of H is 3/4 x 3/4 or 9/16, and the probability 
of H is 1 - 9/16 or 7/16. If C is true/the probability 
of H is 1 and the probability of H is,0. The probabilities 
of the four compound events / found by multiplication, are 
given at the bottom of the diagram. 



9/16 



1/4 




H H 



7/36 5/9 



Figure 3. Diagram for ExatnpLe I. 

The probability of event H, that Mr. and Mrs. A's two 
children are healthy, is 1/4 + 5/9'or 29/36. 

^In more conventional, notat ion, 

P(H) = P(C)'P(H|C) + P(C)-P(H|CK 

^ = 4/9*9/16 + 5/9*l'= 29/36. 

Example 2. Although Mr, P ! s mother died of Hunting- 
ton's disease, Mr. P is now'forty' years of age and healthy. 
If three-qua*ters of all carriers of Huntington's disease 
develop syjrijSforas by the time j;hey are forty, what is the 
probability\hat Mr. P will eventually develop- Huntinton ' s 
disease? 

This is another problem in conditional probability, 
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because there are two conditions under which the son of 
a woman wKl* died of Huntington's diseas'e can be healthy 
at forty. Hither he did not inherit the' abf&rmal allele, 
or he did inherit the abnormal allele but has not yet dev- 
eloped the disease. Both of these possibilities must be 
considered. "S^ 

- In Figure 4,"c represents the event that "Mr. P is 
a carrier, C that he. is not. The probability for each of 
these is 1/2. H is the event that Mr. P is healthy (that 
is, .free of Huntington's disease) at forty, H that he is 
not,. Since three- fourths of all carriers develop the 
disease by age forty, the probabilities of these in the 
event of C are 1/4 and 3/4 respectively. Fn the event of 
C, the respective probabilities are 1 and 0. 




1 

«. Figure 4. Diagram for Example 2. 

Since 1/8 + 1/2 or 5/8 of the children of a woman 
wno died of Huntington's disease can be expected to .Be 
healthy at forty, and 1/8 to be healthy at forty and car- 
riers, the probability that, Mr. P is a carrier given that 

1/8 • 
he is healthy at forty is or 1/5. In 'conventional no- 
tation , 



P(C H) 



. P(H,C) _ 1/8 a 
" P(H) 578 



1/5, 



Exercise 5 . .Bill's father »% mother had Huntington's disease. Although 
half the carriers of Huntington 1 s disease show symptoms by age thirty- 
fitfe, Bill's father just celebrated his thirty- fifth** birthday and is 
still healthy. If Bill has a child someday, what is the probability 



*. that this child will get Huntington's disease? 

Exercise 6 . Willy's maternal grandmother has Huntington^ disease. 
His mother is thirty-six and in good health. If 60% of the carriers 
of Huntington's disease develop symptoms by the time they are thirty- 
• six, what are the chances that Willy wiirt>e free of the disease? 

Exercise 7 . Although both Mr. and Mrs. D have siblings who are Al- 
binos, all three of their children are normal. What are the chances 
that their fourth child will be .an albino? 



4. THE INCIDENCE OF HEREDITARY DISEASE 

4.1 Pat£ern£*<xf Inheritance 

There are now 2000 recognized genetic diseases. 
If tltey are to be investigated and controlled, it is essen- 
*H'aK^ha,t doctors arid^counselors understand the mechanisms 
by' which they a>e inherited. In some cases these mechanisms 
are very ccA^l ex,- involving the interactions of large . 
num k£r s .of jaifferent genes. Dominant and recessive dis- 
orders, hQwever, exhibit simple patterns of inheritance 
that make them easy to identify. 

I 

L In the case of a dominant disorder / every affected 
individual wili have an affected parent. Further, as^'we 
have seen, the' chances that the child of an affected,, parent 
will inherit the, disease are one in two, assuming that* the 
other paren,tq||. normal . ^ In a large group of children of 
such marriages, therefore, close to 501 will get the dis- 
ease^. 

4 -2 Predicting the Incidence of Hereditary Disease with 
Binomial Probability ' ' — ~ 

What can be said about a single family? Suppose", 
that a normal individual marries a carrier of Huntingtqn ! s 
disease and they have three children. What are the pos-" 
sible outcomes? , The diagram in Figure 5 can^ answer^this 
question. * Here K{H) indicates the birth of a normal child 
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with a probability .of 1/.2 ajid A(*s) indicates the birth of 
an affected child with this same probability. To get* the 
probabilities in tHe right*hand column, the individual 
probabilities have been multiplied together," 

Of- the eight possible outcomes, one is for three 
* i * 

* normal- children with a probability of 1/8, 'three, are for 

two normal children and one .affected child with V combined 
probability of 3. * 1/8 or 3/8, three are for one normal 
child and ,two affected children, also with a combined pro- 
bability of 3/8, and one is for'three. affected children 
with a probability' of 1/8, Note that these same results 
could more easily have been obtained by .multiplication 
since 

k (bN> VV) 3 = C(3,0).OsN) 3 (JsA)° + C (3 , 1 ) ftN) 2 (%AJ + 
C(3;2)(^N)(HA) 2 + C(3,3)(^N)°(^A) 3 = 



1/8N J 



3/8N 2 A + 3/8NA 2 + 



1/8A* 



Here the coefficients of the product give the probabilities 
of the corresponding combinations of N ! s and A's. Because 
the binomial theorem is* used to. oftain these results, the 



problem is said to involve binomial probability. 



first child 



N(« 



AOs) 




NNN(l/8) 
NNA(l/8) 
NAN(l/8) 
NAA(l/8) 

ANN(l/8) 
ANA(l/8) 
AAN(l/8^) 
AAA(l/8) 



Figure^. Possible outcomes for a' three child family in 
^ - which one Parent carries- a dominant disease. 

For recessive diseases the probabilities differ.. A 
child will be afflicte'd with a recessive disorder only if 

* " 10 

40 

I. ' ' \~\ 



both of his parents are carriers. When this is the case, 
he has- one chance in four of inheriting* the disease and 
three chances in four of being normal . The possible . 
outcomes for the three child family of two carriers o ( f 
a recessive disorder are shown in Figure -6. 



first child 



H(3/4) 



second child 
.N(3/4)- 



A(l/4) 



Figure 6. 




third child 

N(3/4) 
A(l/4) 
N(3/4) 
A(l/4) 

N(3/4)' 
A(l/4) 
N(3/4) 
A{l/4) 



NNN(27/64) 
NNA(9/64) 
NAN(9/64) 
NAA(3/64) 

ANN(9/64) 
ANAC3/64) 
AAN£3/64) 
AAA (.1/64) 



Possible outcomes for a three child family in 
which both parents carry a recessive disease. 



Of 'the eight possible outcomes, one is for three nor- 

< 

mal children with a probability of 27/64, three are for 
two normal children and one affected child with a combined 
probability of 3 * 9/64 or, 27/64 , three are for one nor- , 
mal child and two affected children with a»combined pro- 
bability of 3 * 3/64 or 9/64, and one is for three affec- 
ted children with a probability of 1/64. As before, a bi- 
nomial expansion would have produced these same results, 
in this case the expans4on of (3/4N + 1/4A) 3 . 

By using the binomial expansion in this ( ^ay', a variety 
of questions can be answered about the incidence of heredi- 
tary disorders. Consider a large ^number of families of 
*four children in fohich both parents^are carriers of cystic 
fibrosis. In what percent of these families will two or-* 
more .'children be affected? 

Solution : The expansion of (3/4N + 1/4A) 4 yields 

C(4,0)(3/4N) 4 (1/4A)° + C (4 , 1) (3/4N) 3 (1/4A) 

+ C(4,2)(3/4N) 2 (1/4A) 2 * 
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+ C(4,3)(3/4N)(1/4A) 3 
+ C(4,4) (3/4N)°(l/4A) 4 . 

» • i 

C 

These last three terms correspond to two, three and four 

affected children respectively. Summing their coefficients 

the probability of two or more children being affected is 

found to be 
a 

C(4,2) (3/4) 2 (l/4) 2 + C(4/3)(3/4)(l/4) 3 + C(4,4)(l/4) 4 

54 + 12+1 _ 67 
\ 256 " 256 ~ ' 

so two or Snore children will be affected in about 261 of 
the families 



Exercise 8 . A man with Huntington's disease has four ; children, 
a) What is the probability that all four will develop Huntington J s 
.disease? b) What is the probability that two children will be af- 
fected and two normal? c) That no more than one will be affected? 

Exercise 9 . , Two parents are carriers of the same recessive gene for 
deafmutism. a) What is the probability that their first child will 
be- deaf ? b) That all five of their children will be normal? c) That 
two of their five children will be^ affected? d) That at least two will 
be affected? 

Exercise 10 > A statistical study is done on families of five children 
in which both parents are carriers of the same recessive disorder. 
In what percent of these families should no more than two children 
be affected? 

Exercise jpf T^ /TcoEple's first child!has cystic fibrosis, a) What 
is the probability that their next child will have it too? b) What 
is the probability that of their next 9 three children only one will 
be affected? j 
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5 . SAMPLE EXAM 



The parents of both Mr. J/iind Mrs. L were carriers of cystic 
fibrosis. What is the probability that neither Mr. L nor Mrs. 
L is a carrier? That both of them are carriers? That both of 
them are carriers but they will have two normal children? 



90% o£*4k£ ^carriers of -Huntington ' s disease haVe symptoms by 
age forty-five. If a man whose father died of Huntington's dis- 
ease is healthy at forty- five, what is the probability that 
he is not a carried 

A study is made of three groups of families in which both parents 

are carriers of a recessive disorder. The families in the first 

group have two children, those in the second group have three, 

and those, in the third have four. If the three groups are of 

equal s£ze, in what percent of the total number of families 

should two or more children be affected? 

i 



i 
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6, ANSife TO EXERCISES 



P(each parent is a carrier) = 2/3. 

P(the child of two carriers is affected) = 1/4. 

2/3 x 2/3 x 1/4 = 1/9. 

P(first child will be bom with the disease) = 1/4. 
P(first two-children will be healthy) = 3/4 x 3/4 = 9/16. 
P(first three children will be carriers) = 1/2 x 1/2 x 1/2 = 1/8. 

P(father is a carrier) = 2/3 so P(first child is an albino) = 

2/3 x 1/2 = 1/3. If first child is an albino, father is a carrier 

and P(each .subsequent child is an albino) = 1/2. 

P(grandfather is a carrier) = 2/3. Assuming his spouse is of 
genotype AA, P(each of hiV children is a carrier) = 1/2 and 
P(each of his grandchildren is a carrier) = 1/4. P(two carriers 
produce a defective' child) = 1/4. 
2/3 x 1/4 x a 1/4 x 1/4 = 1/96. 



P(Bill's father is a carrier) = 



1/4 



= 1/3. 



1/2 + 1/4 

P(Bill is a carrier) = 1/3 * 1/2 = 1/6* - 

P(Bill's child will be a carrier) = 1/6 x 1/2 = 1/12. 




© 1/5 



H H 
3/10 1/2 

P(Willy»s mother is a carrier) = + 1/? 

P(Willy is a carrier) .= 2/7 x 1/2 = 1/7. 
P(Willy is not a carrier) = 1 » 1/7 = 6/7T, 



1/5 



2/7. 
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P(Mr. and Mrs. D are carriers given that they have three normal 
..... . 3/16 27 

children) * 3/16 i 5/9 ; --jjy- / ' • 

P(fourth, child will be ata albino) = 27/107 x 1/4 = 27/4-28 or 
approximately .06.. 

a) C(4,4)ft) 4 '= 1/16 

b) C(4„2)ft) 4 = 6/16 = 3/8 

c) P(one will be affected) = C(4,^%) 4 = 4/16 " 
P(n0^will be affected) = C(4,0)(%) 4 = 1/16 * 

P/no more than one will he affected) s 1/16 ♦ 4/16 = 5/16 
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b) C (5,0) (3/4 r = 243/1024 

c) C(5 t 2)(3/4) 3 (l/4) 2 = 270/1024 



C(5,l)(3/4) (1/4) = 405/1024 
C(5,0)(3/4) 5 = 243/1024 



d) P(none wi»ll be affected) 
P(none will be affected) 
R(at least two will be affected) = 1 - (405/1024 +• ^43/1024) 

= 376/1024 

P(two will be affected) = C(5, 2) (3/4) 3 (1/4) 2 = 270/1024 * 

P(one will be affected) = C(5,l) (3/4) 4 (l/4) = 405/1024 

P(none will be affected) = C(5,0)(3/4) 5 = 243/4024 

270 + 405 +C%3 j 918 _ . 
1024 ~ 1024 " " 



11. a) 1/4 



b) C(3,l)(3/4r(l/4) = 27/64 
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7. ANSWERS TO SAMPLE EXAM 



P(neither Mr. L nor Mrs. L is a carrier) = 1/3'x 1/3 = 1/9. 
P(both are carriers) = 2/5 *. 2/3 = 4/9. 
P(both are carriers and they have two jrtormal children) 
= 4/9 x 3/4 x 3/4 = 1/4. 




H H H H 

1/20 9/20 ) 1/2 0 

PJhe is not a carrier given tpat he is healthy at forty- five) 



1/2 + 1/20 
,2 



10/11 



cte,2)(l/4) = 1/16 " y 

C(3-,2^(3/4)(l/4) 2 * C(3,3)'(l/4) 5 = 9/64 + T/64 = 10/64 

C(4,2f(3/4) 2 (l/4) 2 ♦ C(4,3)(3/4K1/J) ? + C(4,4)(l/4) 4 

<J " 547256 ♦ 12/256 ♦ 1/256 =^ 

* ( J^- 256 * 

1/3 * 1/16 + 1/3 x 10/64 + 1/3 x 67/256 = 123/256 = .48 
.so two or more children should b3 affected in about. 48% o*f the 

& & « 6 

families. 
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BACKGROUND INFORMATION^ 



The bee'js honeycomb is one of the most beautiiul 
geometric configurations in natdSfe. Its Indicate 
geometric form and structural ^onomy have been objects 
of scientific investigation for , centuries . 
* * f. 

The architecture of the . ,r ideal H Tiorieycomb has the 



following basic features: y I 

a) It consists oi: two 'stacks of congruent tubes, 

b) Each tube -is open at .one end and closed at 
the other* 

c) The .two stackjs^of tubes are joined together at 
their closed ends. 

d) Each tubfe has an hexagonal cross 'section . 

e) The closed end of each tube consists of three 
end faces that are congruent equilateral 
paralellograms joined as in Figure 1. 



K Tube Center Line 

End Face ' 




Figure 1. The plane surfaces and surface angles of an ideal 
honeycomb tube*. 
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>Th\ measures of the interior, angles of the 
parallelogf am face's are approximately H0 o 
SsH **^nd 70°. These same angles occur on the *"* 
trapezoidal faces of *tne tubes. Finally, the 
acute angle between the center line of the tube 
and the en'd^/faces is approximately 54°. 

Of cours'e, the^ honeycombs actually constructed**©^ 
wax by the bees deviate from the design describe/^DOve 
in a variety^of ways; The tubes are not exactly %he£$ 
same size arid^ the edges and faces are often curved a 
bit. However, these deviations from the basic plan 
are remarkably small so that the design features listed 
above conform quite closely to reality. 

Why and how do bees follow this exacting design 
pattern for the construction of their honeycombs? As 
you might expect, these are difficult questions to answer 
Some observers have attributed the remarkable uniformity 
of honeycombs to Divine Guidance, others to gerfetic 
natural selection and still others to the action of 
physical forces . 

There' is 'strong evidence to* ,v support the conclusion 
that the architecture of a honeycomb is not merely a 
random selection from a; list of possible multiple tube 
structures of this sort. This evidence is based 
largely on a number 'of theoretical properties of the 
honeycomb structures described above that suggest that 
it is the f, best possible 11 choice -from several points of* 
view. In this paper, we wiijL 7 discuss one of these 
theoretical properties in detail and mention^ some 

others briefly, , ^ . * * * 

* - ■* 

As we pointed out above, ttfe acute angle between 
the center line of tfte^tube in a honeycomb and each of 
its end faces i& approximately 54°, , Why 54°? To provide 
at'least one plausible reason for,: thi*s choice of angle, 
suppose that We assume that the bees construct thefeubes ' 

: • . ' • *g.i 2 
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so that as little wax as possible is used to store a 
given amount of honey.. We could then seek the angle 
between the tube center ]&DJ&*ind end faces tfiat would 
, achieve this minimum.. More precisely, we could pose 
the following theoretical problem: Suppose that a 
hexagonal tube is capped at one end with three equilat- 
eral parallelogram fades and that the other end i§ open. 
Find the angle between the tube center line and end 
faces that will minimize the surface area of the tube 
for a given volume. We wijll solve this problem below 
# and we will see that the correct theor^ical value 
for the angle is approximately 54.7°. This remarkable 
-^agreement with the approximate measured angle of 54° 
^Strongly suggests that the bees are "doing the right 
thing" by some mechanism, or other! 

Scientific observation and stud>^ of the geometric 
structure of the honeycomb can be traced back at least 
to Pappu,s of Alexandria near tfie end of -the third century 
A.D., In thp seventeenth century A. D. , Johannes* Kepler 
was able to explain* the hexagonal cross-section and end 
face structure of theMioneycomb tubes on the basis of 
'packing tubes sub'ject to internal pressures. (This 
explanation will be outlined in the remarks in Section 
3 of this paper.) Jn 1712, Miraldi made some accurate 
measurements df the angles on the caps of the honeycomb 
tubes and conjecture^ that the measuredVa ? ngle of approx- 
imately 54° was the Solution to the minimum problem 
stated in the l^st paragraph. This conjecture was 
verified shortly' later by Koenig. 

2. THE MATHEMATICAL ANALYSIS OF THE BEE'S CELL CAP 

For the purposes of our analysis, we shall regard j 
the individual tubes in a honeycomb to be nexagonal 
cylinders, open at one end and capped at t^e other end 
with^three parallelograms with equal sides as in Figure 2. 
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Figure 2, The tube is a regular hexagonal cylinder, open at 
one end, capped at the other end by congruent 
parallelqgraras. 



\ 



This solid figur-fe can be constructed in the following* 
way : r~ < - .* 

1„ Begin with a solid hexagonai cylinder with end 
faces perpendicular to the center line of the 
cylinder. (See FigurG*"^-! 




Figure 3. A solid, regular, hexagonal right cylinder. 



2) Let A and -B. denote the upper left-hand and upper 
right-hand corners of two adjacent lateral faces 
and let C be-a point down the edge joining these 
faces . 

3) Pass a plane through Uhe points A, B, C to cut 
off a tetrahedral "corner 11 ABCD of hexagonal 
cylinder. ' 

4) Using the line segment, joining A and B as a 
hinge, swing the tetrahedron upward through 
anangle of 180^. (See Figure 4.) 

i C 




:Figure 4. Tetrahedron ABCD relocated as tetrahedron ABCD'. 

Jhe point C will move to the new position X 1 
and D to a new position D'^ The- points A,C,B,C ! 
will be successive vertices of a parallelogram 
and C, D* will lie on the axis of "the cylinder. 
(These statements can be verified by checking > 
the angles involved, in the tetrahedral corner. 
You should carry out the details!) 
5) j Repeat the above procedure' with the remaining 

two pairs of adjacent faces to obtain the solid / 
in Figure 1. 9 

Suppose that 0 is/ the acute angle at C 1 between the 
center line* of the hexagonal cylinder and the face ACBC 1 . 
(By construction, 8 is also the acute angle at C between 

the verticalv^dge of the hepc^gonal cylinder and the .face 

■„*> 



ACBC.) Let h dentfte the height of the^rlginal cylinder 
(before the three tetrahedral corners were cut) and let 
s be the*width of each face. 




Figure 5. The dimensions of the cylinder as it was before the 
tetrahedral corners were cut off and moved. 

Then, by virtue of the construction described above, 

the volume of the hexago^a^tyl inder capped with the 

three equilateral paral leldgram ^faces is equal to the 

volume cylinder of height h and side width s. Since thej* 

^entral angle subtended by. each face is tt/3, it follows 

that . • <0 

/o . , 
7 s]h 
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Next, we shall derive an expression for the total 
surface area of the capped hexagonal cylinder. The 
diagonals M and CC 1 " of the parallelogram face ACBC' / 
bisect each other at a point E and, since ACBC 1 has 
e^ual sides, at fallows that ACBC 1 is divided into four 
c6ngruent right triangles with side length's equal to 
the lengths of AE and CE. 



One end face - 

-Tube Center Line 




Figure 6. The quadri latest AC BC 1 is a parallelogram whtase 
diagonals meet at right angles. 

Notice that the length of AE is s • (/J/2) and that the 
length of CE is ts/2) • esc 9 . Consequently, the area of 
ea.cn face of the ; cap is ' * • 



/S CSC 8^ 2 [JZ a 

( 2 ) s [ T^ 0 ; 6 



and so the total surface area of the cap is 

(2) ^(s 2 csc e ).■ . 

The capped hexagonal cylinder has six congruent trape- 
zoidal faces. The width of these faces is s, one vertical 

edge has height h and the other has height 

/ 

h - length UP = h - |cot & 

Therefore, the area of each lJteral face) is 

2 ' 

s • ^(h + [h - -|cote]) = hs*- -coto 
and so the < total surface area of the lateral faces is 

' ' 3s 2 • ' 

(3) 6hs- - ^-cot 6 • 



We now add (2) and (3), to obtain the total surface area 

S of the capped hexagonal cylinder (open at the other end^ : 

»■ • '. 

7 




53 



(4) S = 6hs'+ |s 2 (-cot 6 + /3 esc 6 ) 



Thus, the surface area S depends on three variables 
h, s, 9. If, for a fixed' cfioice^of n and s, we seek to 
minimize S as a function of 6», we differentiate (4) with 
respect to 6 .and equate the result to (t \6 *f ind the 
critical points : % » 

(5) |s 2 (csc?.a^- /IcscedotG) = 0 

(6) esc 6 - /5cot e = 0 *„ 

l7) cos ; =£21a = ^l: 

J esc e > *j • 

* * t • 

The acute angle satisfying th.e last equation is 1 -) 

f«) e 0 = cos' 1 ^) : 54.70. 

The First Derivative Test can be used to check that this, 
critical point is actually a minimum. **j k ' 

It is very important to notice at this poin£ that 
the value of 8 that yd e lds tn ^ minimum surface area S 
does not depend on the variables h and s. fn other 
words, although w£ fixed s and h before we minimized 
S with resplct to (J the value of'Q of approximately 
34. 7° that minimized the surface area S is the same for 
zny oh^-ige of the dimensions a and h. This is actually 
a stronger conclusion than the one we set out to „ 
establish: For a given* volume, the value of 8 that 
minimized the surface area S of the hexagona I cy Under 
zapped at one end with three equilateral parallelograms 
is approximately 54,7°, 

3. SOME ADDITIONAL COMMENTS 



Some of the other features of the geometric 
structure of honeycombs can be explained on the^basis 
of internal forces within the tubes. We shall now 
outline these ideas briefly. * • 

CO 



.8 




Suppose we are given a system of circular cylinders 
that is closely packed so that each cylinder is in con- 
tact with six others along lines parallel to the center 
line of the cylinders. A cross : sectional view of a 
central cylinder and its immediate neighbors is given 
in Figure 7. 




Figure 7. Cross section of seven congruent circular cylinders 
packed together. * 

- - ■ . ' ^ 

Suppose that the system is subjected to a uniform 
outward pressure from within each cylinder. This 
pressure will push th^ cylindrical walls out Into the 
empty spaces between the cylinders. It can be shown 
'that. when all this empty space has been filled, the 
system will consist of hexagonal cylinders.' Thus, it 
j.s quite conceivable that the hexagonal cross-sectional 
structure of the honeycomb is due to the efforts of the 
bees to pack as much honey as possible in each tube. 

Jt A similar explanation can be given^for the geometric 
Sffucture of the caps of the tubes in a honeycomb. Consider 
a system of two stacks of circular cylinders, each 
cylinder open on one end and capped on the other with a 
hemispherical cap like attest tube. Plac^the stacks so \ 
'* '+ ' 




si 



that the capped endsVof the cylinders in, one stack are 
in close contact with the capped ends of the other stack 
In the closest possible arrangement of the two stacks, 
each hemispherical end of a cylinder in one stack is 
in contact with three hemispherical ends of cylinders 
in the other stack. 

If /his system is subjected to uni fornix outwarcf 
pressur4 from within each cylinder, the walls of the 
cylinders push out to fill t|e empty space between the 
cylinders. We have already pointed out that the final 
lateral surfaces of the cylinders will have equal 
hexagonal cross - sections . It can also be shown that 
the end^surfaces will consist of three equilateral * 
parallelogram faces when the empty spaces at^the ends 
are completely filled. Thus, the .basic geometry 
configuration of the honeycomb could result from the 
action of a uniform internal' pressure. 

Do the bees construct the j^st possible caps for 
their honeycomb tubes? We have shown that they do'if 
the caps are required to consist of three equilateral 
parallelogram faces'. But are there "better choices of 
polyhedral surfaces for the cap? The answer is: Yes! 
for example,, it has been shown that a cap consisting 
of two hexagonal and two equilateral parallelogram 
surfaces joined as indicated in Figured: 




V Figure 8. # Honeycombs shaped like this would be slightly more 
- efficient. " ' 
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would be mote efficient than a cap with three equilateral 

parallelogram* faces . However, the reduction im surface 

area per unit volume for this more complicated cap is 

very small. The bees probably decided that these fancy 

caps are not cost effective when both labor and materials 

♦ 

are considered! 

* « 

4. PROBLEMS 

1. Verify the^foTi owing, details of the construction described \x\ 
figure 4: W ^ 

a) *nje points A, C, B, C are successive vertices of a plane ^ 
paralellogram. 

b) C, D' lie on the axis of the hexagonal/cylinder . 

2. Carry out the details of the First Derivative Test to verify ' '= 

irtiat *L = 54.7° is actually a minimum of -the surface area 
\ J 0 

Mhct ion. (See Equations (5), (6), (7).) 

3. Among all regular polygons, only the square, the equilateral 
triangle and the hexagon can be used to "tile" the plane 
(that is, to cover the plane wjtth congruent non-overlapping y 
pieces)\ Show that of these three polygons, the hexagon has w 
the least perimeter for a given enclosed area. 

5, SOME REFERENCES FOR FURTHER READING 

1. Thompson, D'Arcy W. , On Growth and Form Cambridge 
University .Press , 1917. -(Available in the Mat hematics^ 
Library under the call number 577.3 T 37o.) This book 
is "a very readable and fascinating treatise on geome- 
tric forms in nature. 

2, T6-th, Lv Fejes, What the bees know and whaj^ they -do ^ 
not know; Bulletin of the American Mathematical 
Society, Vol 70 (1964) pp 468^481. (Mathematical 

* Library call number: 510,6 Am B2,) This paper is 
the text of an invited address presented at a national 
meeting of the ^American Mathematical Society in 196^.- 
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It is a very readable presentation of some minimum surface 
area problems ^Yel ated to honeycombs. In particular, a 
full discussion of the comparative economy of the tube 
caps in Figures- 1 and 8 is given. 

' 6. SOLUTIONS TO PROBLEMS 

t 

1. (a) Since the point D is rotated through an angle of 180° to 
the new position D', the points A, B, C and C lie in a plane. 
The triangles ABC and ABC are congruent by construction, so 
A, C, B, C are- successive vertices of a plane parallelogram. 

1. (b) Sin<>e angle )ADB = 120°,,\it follows that angle }AD'B = 120° 
so D ' is on the axis of the hexagonal cylinder: , Since the line 
segment DC is parallel to this axis and since the segment DC is 
rotated through 180° to the new position D*C, it follows that 
D'C is parallel to this axis also. Since the point D 1 is on 
this axis,- so is the point C , 

2. Note that 

* * 
ds 3 2 ' - * ' 

d6 = 2 S ( csc " ® " ^3 csc 9 _ cot 6) * 
= j s 2 cs^'e (esc 0 - cote). 

' 3 D * Tf 

Since j s^>0 and esc 9 > 0 for 0 < 6 < y , it is only necessary to 
show that "esc 6 - ^3 cot 8 changes from- negative to positive 
vafues as 6 increases through 8 Q = cos _1 (l//3). This sign 
change can be verified by inspecting the graphs of esc 8 and 
/3 cotG., 

3. Suppose that the area enclosed by the polygon is A^ that the 
side length is s, and that the perimeter is P. - Then for the: 

(a) square, A = s 2 , P = 4s \ 

(b) equilateral triangle, A = — j — > P = t 3s 

3/3 s 2 * * 

(c) regular hexagon, A = — ^ — , P = 6s 

If we express P in terms of A, we obtain for the 
£*X square, P = 4^ 
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(b) equilateral triangle, P = /K ^ 
•(c) regular hexagon, P = a -a /K 

* * 

From these calculations r *oae can conclude that for a given 
area A, the perimeter of the equilateral triangle is greater 
than that of the square which in turn is greater than that of, 
the regular hexagon » • . # 
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f STUDENT FORM 1 

Request for Help 



Return to: \ 

EDC/Ul^ 

55 Chapel St*V 

Newton, MA 0iM60 



Student : If- you have trouble with a specific part of this unit, please fill 
6ut this form and take it to your instructor for assistance. The information 
you give will help the autho? to revise the unit. 



Your N&me 



Unit No. 



Page_ 



GpUpper 
OMiddle 
,0 Lower* 



OR 



Section 



Paragraph^ 



OR 



Model Exam 
Problem No._ 

Text 
Problem No. 



Description* of Difficulty: (Please be specific) 



Instructor : Please indicate your resolution of the difficulty in this box. 
Corrected errors in materials. List corrections here: 



o 



Gave student better explanation, example, or procedure than in unit. 
Give brief outline of your addition here: 



o 



Assisted student in acquiring general learning and problem-solving 
skills (not using examples from this unit.) 



es 



Instructor's Signature^ 



Ple^aSe^use reverse if necessary. 



Return to: 

STUDENT FORM 2 EDC/UMAP 

Unit Questionnaire , , £ 5 ^ pe l A S !;- 

x Newton, MA 02160 



Name ^ Unj.t No. Date 

> I ~ — 4 ; — 

Institution * » Course No. I 



Check the choice for each question that comes closest to your personal opinion. 

1. How useful was the amount of detail in the unit ? N 

Not enough detail to understand the unit 

Unif would have been* clearer with more detail 



Appropriate amount of detail 

Unit was occasionally too detailed, but this was not distracting 

< Too much detail; I was. often distracted 



2 . How helpful were the problem answers ' 

Sample solutions were too brief; I could not do the intermediate^teps 

Sufficient inf oljflfeion <was given to solve the problems 

lS^Br€ 



Sample solutions^Bre too detailed; I didn't need them 



3<* Except for fulfilling jthe prerequisites, how much did you use other sources (for 
* example, instructor, friends, or other books) in order to understand the unit? 

A Lot Somewhat A Little Not at all 



4. How Xong was this unit in comparison to the amount of time you generally spend on 
a lesson (lecture and homework assignment) in a typical math or science course? 

'Much Somewhat About Somewhat Much 
_Longer , Longer the Same 'Shorter Shorter 

— 7 

5. Were any of the following parts of the unit confusing or distractjdlg ? (Check 
as many as apply.) ^ 

Prerequisites 

Statement of skills and concepts (objectives) f 



Paragraph headings 
^Examples j 
^Special 'Assistance Supplement (if present) 

Other, please explain I 



6. Were any of the following parts of the unit particularly helpful? (Check as many 
as apply . ) « 

Prerequisites ^ 

Statement of skills and concepts (objectives) . % 

Examples 



Problems 
— * 

Paragraph headings 



_Table of Contents 

^Special Assistance* Supplement (if present) 
Other, please explain 



Please describe .anything in the unit that you did not particularly like. 



Please describe anything that you found particularly helpful. (Please use the back of 
this sheet if you need more space.) _ I 
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